Density scaling in viscous liquids: Prom relaxation times to four-point susceptibilities 
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We present numerical calculations of a four-point dynamic susceptibility, X4(^)i for the Kob- 
Andersen Lennard- Jones mixture as a function of temperature T and density p. Over a relevant 
range of T and p, the full t-dependence of X4(^) and thus the maximum in X4(i)) which is proportional 
to the dynamic correlation volume, are invariant for state points for which the scaling variable p^ jT 
is constant. The value of the material constant 7 is the same as that which superposes the relaxation 
time, r, of the system versus p'' /T. Thus, the dynamic correlation volume is a unique function of 
r for any thermodynamic condition in the regime where density scaling holds. Finally, we examine 
the conditions under which the density scaling properties are related to the existence of strong 
correlations between pressure and energy fluctuations. 
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During the last decade extensive evidence has accu- 
mulated that the dynamics of molecules in supercooled, 
highly viscous liquids is inherently heterogeneous 

' that is, spatial variations in mobility persist for times 
commensurate with the structural relaxation time r. 
Thus, there is growing interest in characterizing the dy- 
namic heterogeneities to better understand the phenom- 
ena associated with the many-body dynamics of vit- 
ryifing liquids. These phenomena include rotational- 
translational decoupling [1, Q, the dynamic cross-over 

! @i 13' enhanced mobility under confinement @, non- 
exponentiality of the relaxation [9!], and most promi- 
nently, the slowing down of the translational and rota- 
tional dynamics upon approach to the glass transition 
temperature Tg. To explain these features, theories of the 
glass transition as diverse as the classical Adam-Gibbs 
entropy [iO|] and Cohen- Grest free volume [ll| models, 

' as well as more modern approaches [l^ IT3| , invoke dy- 
namic heterogeneity having a length scale that grows in 
concert with r upon drawing near to Tg. 

A proper description of the dynamic heterogeneities re- 
quires multi-point dynamic susceptibilites, which reflect 
correlations in the spatial variation of the dynamics. A 
four-point dynamic susceptibility Xi{t) can be calculated 

' as the variance of the self-intermediate scattering func- 
tion F,(/c,i) 

I 

Xi{t)=N[{f,(k,t))-F^{k,t)\ (1) 

where fs{k,t) is the instantaneous value, such that 
{fs{k, t)) — Fs{k, t). X4:{t) quantifies the amplitude of the 
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fluctuations associated with fs{k, t) and has a maximum, 
X™"^ — X4(trr,.an-), proportional to the dynamic correla- 
tion volume [3ll3- X4(0 exhibits various regimes [3], 
but most interesting for study of the glass transition is 
the behavior around tmax ~ t [iB] ■ Recent numerical and 
simulation works have shown that the dynamic correla- 
tion volume grows upon cooling [13, [H, [l^. However, 
very little is known about the combined temperature and 
density dependences of X4(^)- particular, unexplored 
is the possibility of a description of the p and T depen- 
dences of X4(0 i'^ terms of the scahng pro per ty estab- 
lished for the structural relaxation time [20, l2l|, I24I 

t = J-i(pVT) (2) 

where Ti is a function and 7 a material constant. Experi- 
ments have shown that Eq. ^ , with similar relations for 
the diffusion constant and viscosit y, a pplies universally 
to organic, non-associated liquids [23|, with a range of 
validity extending from the high temperature Arrhenius 
regime down to Tg [23|. We anticipate a similar relation 
for tmax and examine the possibility that this scaling 
property extends to the dynamic correlation volume 

x'r''=WlT) (3) 

To cast this work in more general terms, we assess 
whether for a prototypical model glass-former, the den- 
sity scaling properties hold for the full time dependence 
of both Fg{k,t) and X4(i), using the same scaling expo- 
nent 7. The connection between the dynamic properties 
of viscous liquids and 7 is intriguing because molecular 
dynamics (MD) simulations have shown this parameter 
to be a measure of the correlation between fluctuations 
in the potential energy and the virial [1^ [2^ [l^] , pos- 
sibly reflecting a hidden scale invariance in viscous liq- 
uids dl]. 7 is also related to the steepness of the effec- 
tive repulsive potential in the range of closest-approach 
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FIG. 1: Reduced relaxation times r* = r^p^^^T^^^) as a func- 
tion of p"'/r with 7 = 5.1 for all studied densities. Inset: 
self intermediate scattering functions as a function of reduced 
time t* = t{p^''^T^'^) for state points at which p^/T = 5.07: 
T = 0.402 (at p = 1.15), T = 0.50 (at p = 1.20), T = 616 (at 
p = 1.30), and T = 0.912 (at p = 1.350). A constant reduced 
wave-vector k* = k(p^^^) = 7.44 is considered. 



between particles [29|, |30|, |3l| , although this connection 
may seem tenuous given the limitations of a two-body po- 
tential in describing interactions in real liquids and the 
possible non-trivial role of attractive forces [sT], [H, [s^ . 
We further analyze this herein by evaluatin g p ressure- 
energy correlations for the model studied in [33], which 
did not conform to density scaling over the relevant den- 
sity regime. 

In this work we consider the well-studied Kob- 
Andersen (KA) binary mixture [s^l as a model glass- 
forming liquid. It consists of 1000 particles in a cubic 
box with periodic boundary conditions. Particles inter- 
act through the Lcnnard- Jones (LJ) potential 



(^)"-(^) 



(4) 



where a, (3 — 1,2 are species indices. The values of the 
parameters in Eq. ^ can be found in the original paper 
[3J|. In the following we use reduced LJ units, assuming 
(Til, en, and miaf^/cn (where mi is mass) as units of 
distance, energy, and time, respectively. We performed 
MD simulations in the NVT ensemble using the Nose- 
Poincare thermostat [s^ with a mass parameter (5=5.0. 
We considered five isochoric paths in the density range 
1.150 < p < 1.350. For each state point we averaged the 
dynamic properties over 20 independent realizations of 
the system. 

Following previous work and to emphasize the con- 
nection with the exact density scaling relations observed 
for inverse power law (IPL) potentials, we analyze "re- 
duced" quantities (indicated by stars) using p~^l^ and 
T^l"^ as reduction parameters for distances and veloc- 
ities, respectively. We begin our investigation of the 



FIG. 2: Maximum of four point dynamic susceptibility as a 
function of p~' jT with 7 = 5.1 for all studied densities. Inset: 
four point dynamic susceptibility as a function of reduced time 
t* for state points at which p^/T = 5.07 (same state points 
as inset of Fig. [1} . 



density scaling properties of the KA model by calcu- 
lating f^,(fc,t) = (l/^)Ef=i(exp{«k-[r,(i)-rj(0)]}) 
for the various densities at a fixed reduced wave-vector 
k* = k{p^^^) = 7.44, which matches the position of the 
first peak in the static structure factor {k = 7.0) for the 
well-studied density p = 1.2. The reduced relaxation 
times T*, defined as the time for Fs(k* ,t*) to decay by a 
factor of e, are shown in Fig. [T] versus the scaling variable 
p'' /T. The material constant 7 — 5.1 ± 0.1 provides the 
optimal collapse of r* onto a single curve; this value is in 
accord with the scaling behavior found previously for the 
diffusion coefficient of this system ^30^ . We note that the 
pressures attained at the lowest temperatures range from 
1 to 20 reduced LJ units, depending on p. Assuming Ar- 
gon units, this corresponds to 0.04 to 0.8 GPa, which is a 
significant and experimentally accessible pressure range. 

It has been demonstrated experimentally that the 
shape of the frequency-dependent linear response func- 
tion is a function of r and thus of p^ jT [36!, [37|, a result 
consistent with the existence of "isomorphic" points in 
liquid state diagrams [11]. We calculated the interme- 
diate scattering functions for each of the five densities 
at the respective temperatures corresponding to a fixed 
value of p'^ jT = 5.07. In the inset of Fig. [1] these cor- 
relation functions are plotted as a function of reduced 
time, and, as found previously over a more limited den- 
sity range [111, Fs{k*,t*) has essentially the same shape 
for state points for which p^ /T is constant. Thus, not 
only do the relaxation times superpose as a function of 
p'^ /T , but the entire t-dependence of the correlation func- 
tions is invariant for isomorphic state points. 

We now examine the correlations in the dynamic fluc- 
tuations to assess explicitly whether density scaling can 
be extended to high-order correlation functions, as envis- 
aged in [ll]. To do this we calculate the four point dy- 
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FIG. 3: Maximum of four point dynamic susceptibility as a 
function of reduced relajcation times r* for all studied densi- 
ties p. Inset: reduced time at which X4(i*) is maximum as a 
function of r*. 



namic susceptibility [Eq. ([T])] associated with the complex 
instantaneous value of fs{k,t) at the same fixed reduced 
wave- vector k* considered above. The use of the com- 
plex self intermediate scattering function (rather than 
the real part used in previous studies) removes the finite, 
long-time limit from Xii^) without altering the general 
features of the correlation function. In Fig. [2] the maxi- 
mum of X4(t) is plotted as a function of p'* /T using the 
aforementioned 7 = 5.1. To provide an estimate of the 
statistical uncertainties, we include errors bars that rep- 
resent two standard deviations on the average over sys- 
tem realizations for selected states. As seen in Fig. [51 
density scaling applies to xT"'^ within the estimated er- 
ror, using the same value for the scaling exponent that 
superposes the relaxation times. In the inset of Fig. [2l 
we show the full i-dependence of X4(t) for state points 
at which p'' /T — 5.07. As is the case for the "average" 
intermediate scattering functions, within the estimated 
error bars the Xi{t) f^-H 011 ^ single curve. Deviations 
from the scaling are observed only for the lowest investi- 
gated density (p = 1.15). 

To make contact with previous numerical and exper- 
imental work on four point dynamic susceptibilities, we 
show in Fig. [3] x™"^ as a function of the reduced re- 
laxation times, T* . Since xT""^ is proportional to the dy- 
namic correlation volume, Fig. [3] confirms the presence of 
a steady but rather mild growth of dynamic correlations 
as the structural relaxation times increase to the point of 
vitrification. It also shows that, at fixed t*, the dynamic 
correlations are invariant to either T or p, as expected 
from the density scaling of both the average dynamics 
and the dynamic correlations for the same value of 7. 
Recent experiments |T9| have shown that a similar result 
holds for the temperature derivatives of the two point dy- 
namic correlation function, xt (t) , which provides a lower 
bound to X4(t) [H, |T3|- This correspondence supports 



FIG. 4: Slope of U-W correlations as function of T at p = 1.2 
(circles) and p = 1.3 (squares) for the KA LJ model (open 
symbols) and its WCA variant (filled symbols). The dashed 
lines indicate the range of the dynamic scaling exponent 7 = 
5.1 ± 0.1. Inset: Pearson correlation coefficient as a function 
of T for p = 1.2. 



the validity of the experimentally accessible XT(t) as an 
approximation to Xi{t)- the inset of Fig. [31 we show 
the T*-dependence of the reduced time associated with 
the maximum in Xii^)- ^s expected, the two quantities 
are essentially equal. 

Recent numerical work psi . [26l . \2T\ has shown that the 
dynamic scaling exponent 7 can be independently esti- 
mated on the basis on the correlation between fluctu- 
ations of two thermodynamic quantities, the potential 
energy, f/, and the virial, W (the configurational part 
of the pressure). These fluctuations, AU = U — (U) 
and AW — W ~ (W) are proportional for particles in- 
teracting with IPL potentials [2^ and have been shown 
to be strongly correlated (Pearson correlation coefficients 
R > 0.9) for various other liquids [26,]. In the latter cases 
the slopes, obtained from linear regression of AW vs. 
AU, are equal within the statistical fluctuations to the 
dynamic scaling exponent Such results support the 
conjecture that hquids display strong U-W correlations 
if and only if they comply with density scaling [2^, the 
inference being that these properties have a common ori- 
gin in the same generalized IPL approximation of the 
interaction potential [2^. 

To further elaborate on these aspects, we show in Fig.[31 
the AW vs. AU slopes, F, obtained for the KA model, 
as a function of T for two densities (p = 1.2 and 1.3). 
r(r, p) is always close to the scaling exponent (within the 
estimated uncertainty), as found previously from simu- 
lations along isobaric paths [27]. Also included in Fig.[4l 
are results from additional simulations carried out using 
the purely repulsive, Weeks-Chandler-Andersen (WCA) 
variant of the KA model [H, S^]. In the WCA model, 
the interaction parameters are unchanged but each of the 
pair potentials, Uaf3{r), is shifted so that the minimum is 
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zero and the potential is truncated at this minimum . 
Very recently, Berthier and Tarjus i33|] concluded from 
simulations of this model that density scaling surpris- 
ingly requires the contribution of the attractive interac- 
tions; in the WCA model, in fact, the scaling of r was 
absent except at very high densities. As seen from Fig. [4] 
and its inset, in the WCA model there are strong U-W 
correlations {R > 0.9); however, the value of F changes 
significantly as the state parameters are varied. That is, 
the local scaling exponents change too much (for reasons 
yet unknown), causing a breakdown of the density scal- 
ing. Nevertheless, for every state point U and W are 
strongly correlated. 

In the light of these results, the conjecture of Pedersen 
et al. [l^ connecting strong U-W correlations and den- 
sity scaling must be partly reformulated. The existence 
of strong U-W correlations is not a sufficient condition 
for Eq. ^ to apply. There is an additional requirement, 
that the slopes of the U-W correlation must be (almost) 
insensitive to variations of the state parameters. This is 
consistent with previous simulation and experimental re- 
sults on local dynamic scaling exponents [4l|, \^ . Under 
these conditions, density scaling applies to a very good 
approximation to all time-dependent properties, includ- 
ing high-order time-dependent correlation functions, as 
conjectured in 'ss']. On the other hand, further work is 
needed to understand under which conditions the slopes 



of U-W correlation are insensitive to variations of p and 
T, and how this is related to the attractive part of the 
interaction potential. 

From simulations over a range of T and p, the dynamic 
correlations in a supercooled LJ mixture arc shown to 
have a spatial extent that depends only on the quantity 
p'^ /T. In consideration of the general behavior of liquids 
obeying density scaling, this means the dynamic correla- 
tion volume is related to other mutually correlated prop- 
erties: the relaxation time, the shape of the relaxation 
function, U-W correlations, and at least approximately 
the isobaric fragility. From experiments it is also found 
that the relaxation time is constant both at the onset 
of non-Arrhenius behavior at high T (i^ and at the dy- 
namic crossover at ca. 1.2Tg [3,14^. Since x™""^ depends 
only on r, these changes in the dynamics also occur at 
a fixed (pressure-independent) correlation volume. Al- 
though examining repulsive exponents other than the LJ 
value of 12 remains for future work, if the equivalence of 
the scaling exponents for x™°^ and r is maintained, we 
expect that the dynamic correlation volume should vary 
among different materials. 
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